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Recent Chandra observations have convincingly shown that the soft X-ray emission from the isolated neutron 
star candidate RX J1856. 5-3754 is best represented by a featureless blackbody spectrum, in apparent contrast 
with the predictions of current neutron star atmospheric models. Moreover, the recently measured star distance 
(« 120 — 140 pc) implies a radiation radius of at most ~ 5 — 6 km, too small for any neutron star equation of state. 
Proposed explanations for such a small radius include a reduced X-ray emitting region (as a heated polar cap), 
or the presence of a more compact object, as a bare quark/strange star. However, both interpretations rely on 
the presumption that the quark star or the cap radiates a pure blackbody spectrum, and no justification for this 
assumption has been presented yet. Here we discuss an alternative possibility. Cool neutron stars (T < 10^ K) 
endowed with a rather high magnetic field {B J> 10^^ G) may suffer a phase transition in the outermost layers. 
As a consequence the neutron star is left bare of the gaseous atmosphere ("naked"). We computed spectra from 
naked neutron stars with a surface Fe composition. Depending on B, we found that the emission in the 0.1-2 keV 
range can be featureless and virtually indistinguishable from a blackbody. Moreover, owing to the reduced surface 
emissivity, the star only radiates ~ 30 — 50% of the blackbody power and this implies that the size of the emitting 
region is larger than for a perfect planckian emitter for the same luminosity. When applied to RX J1856. 5-3754 
our model accounts for the observed X-ray properties and solves the paradox of the small radius: we predict an 
apparent star radius of ~ 10 — 12 km, consistent with equations of state of a neutron star. The optical emission 
of RX J1856. 5-3754 may be explained by the presence a thin gaseous shell on the top of the Fe condensate. 



1. Introduction 

The family of thermally emitting isolated neu- 
tron stars (NSs) includes seven peculiar ob- 
jects serendipitously discovered in ROSAT PSPC 
pointings (see e.g. [1] for a review). These 
sources are characterized by similar properties: a 
blackbody-like, soft spectrum with Tbi, 100 eV; 
low X-ray luminosity, Lx ~ lO'^" — lO'^^ergs^^; 
low column density, Nh ~ 10^*^ cm~^; pulsa- 
tions in the 5-20 s range (detected in four sources 
so far). Due to their proximity and to the fact 
that they are "truly isolated" (i.e. not associ- 
ated with a supernova remnant), they are ex- 
cellent candidates to probe directly the emission 
from the star surface. On the other hand, since 
they are dim and extremely soft, until recently 
high quality spectra were not available. PSPC 
data only provided evidence that a hard tail is 
absent in all cases and that a blackbody gives 
a satisfactory spectral fit. The situation has re- 
cently improved for the two brightest and clos- 



est objects, RX J1856. 5-3754 and RX J0720.4- 
3125, which have been target of deep observations 
with Chandra and XMM-Newton ([2,4,5]). The 
hope was to probe directly the neutron star sur- 
face composition, by detecting some of the spec- 
tral feautures predicted by atmospheric models at 
different chemical composition and/or magnetic 
field strenght. Surprisingly/disappointingly, in 
both sources no spectral features have been de- 
tected. In the case of RX J1856. 5-3754, the re- 
sult is particularly robust and striking: a very 
long (~ 500 ks) Chandra pointing has shown that 
the spectrum is better fitted by a simple black- 
body than by more sophysticated atmospheric 
models. Also, upper limits on the pulsed fraction 
of RX J1856. 5-3754 are extremely tight, < 3% 
([5]).^ The sources is known to exhibit an "op- 
tical excess" : accurate photometry of the optical 
counterpart with combined Very Large Telescope 



few weeks before this meeting, J. Triimper reported at 
the 34*'' COSPAR Scientif Assembly an even more strin- 
gent limit, ~ 1%. This result is now published on [6]. 
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and Hubble Space Telescope data has shown that 
the UV-optical emission from RX J1856. 5-3754 
exceeds the Raylcigh-Jcans tail of the X-ray best- 
fitting blackbody by a factor ~ 6 ([7]). The value 
of the distance {d « 120 - 140 pc; [8,7]), used in 
conjunction with Chandra data, yields a radiation 
radius of only ~ 5-6 km ([5]), too small for any 
known NS equation of state. 

The small apparent radius and the blackbody 
X-ray spectrum led to the intriguing suggestion 
that RX J1856. 5-3754 might host a quark star 
([5,9]). The motivation is that the radiation ra- 
dius is compatible only with equations of state 
(EOSs) involving strange matter and that bare 
quark stars, i.e. those not covered by a layer 
of ordinary matter which act as an atmosphere, 
would presumably emit a pure blackbody spec- 
trum. While a quark star is a conceivable op- 
tion, present observations of RX J1856. 5-3754 do 
not necessarily demand this solution. Other more 
conventional model fits, based on two blackbod- 
ies, can account for both the X-ray and opti- 
cal emission of RX J1856. 5-3754, giving at the 
same time acceptable values for the stellar radius 
([10,7,11]). In this scenario X-rays originate in a 
relatively large {0 ~ 20°) heated polar cap, while 
the rest of the surface is colder (~ 30 eV) and pro- 
duces the optical emission. However, since both 
regions should be covered with an optically thick 
atmosphere, the reason why they should produce 
a featureless spectrum has not been thoroughly 
understood yet. 

In this paper we suggest an alternative possibil- 
ity. As pointed out by [12] (see also [13]) NSs may 
be left without an atmosphere by the onset of a 
phase transition that turns the gaseous outer lay- 
ers into a solid. This happens if the surface tem- 
perature drops below a critical value Tcrit which 
in turn depends on the star magnetic field. The 
source then appears as a "naked" or "bare" neu- 
tron star: we only detect the emission from the 
solid layers. While it appears difficult even for the 
coolest isolated NSs to meet the requirements for 
the onset of a phase transition in light element 
(H, He) surface layers, this may be the case if the 
external layers are dominated by heavy elements 
(such as Fe). The uncertainties on the conditions 
for Fe condensation are currently quite large, but, 
depending on the magnetic field, it is possible 
that the surface temperature of RX J1856. 5-3754 
falls below the critical value. 

If condensation occurs, the emitted spectrum is 



not necessarily a blackbody: an overall reduction 

of the surface emissivity and strong absorption 
features are in fact expected when the photon en- 
ergy becomes lower than the plasma frequency of 
the medium. However these absorption features 
may or may not appear at soft X-ray energies 
depending on the model parameters. In the fol- 
lowing we investigate the emission properties of 
naked NSs and discuss the relevance of our model 
in connection with the NS candidate RX J1856.5- 
3754. For all details we refer to [14]. 

After this work has been completed we become 
aware that similar conclusions have been indepen- 
dently reached by [6] and have been reported by J. 
Triimpcr at the 34"^ COSPAR Scientif Assembly. 
These authors attempted a fit to combined UV- 
optical and X-ray data of RX J1856. 5-3754 with a 
depressed blackbody plus an enhanced Rayleigh- 
Jeans tail, by miming phenomenologically the sit- 
uation here described theoretically. Their results 
for the star radius are close to ours. 

2. The Model 

2.1. Bare Neutron Stars 

Theoretical research on matter in superstrong 
fields started over 40 years ago and, although 
many uncertainties still remain, much progress 
has been made especially for H and He com- 
positions (see [13] for a recent review and ref- 
erences therein). At B ^ Bo = rriee^c/fC' ~ 
2.35 X 10^ G strong magnetic confinement acts on 
electrons, and atoms attain a cylindrical shape. 
It is possible for these elongated atoms to form 
molecular chains by covalent bonding along the 
field direction. Interactions between the linear 
chains can then lead to the formation of three- 
dimensional condensates. The critical tempera- 
ture below which phase separation between con- 
densed H and vapor occurs is given by [12] 



crit 



0.1 



IM.lB^f 



4.4 (lnSi2 -6.05)^ - 



eV, 



( 



where B12 = B/{W^^G), the 
proton cyclotron and plasma frequencies. For 
heavier elements (such as Fe) all estimates are 
still quite crude. T^ru is obtained by equating 
the ion density of the condensed phase near zero 
pressure to the vapor density, but all available 
models are approximate near zero pressure. The 
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most recent estimate of Tcrit for Fe phase sep- 
aration, as given by [13], is Tj^fj ?a 21B\^2 ■ 
The density of the condensate (at zero pressure) 
is « 560Z-3/5Bi^^gcm-3 (here Z denotes 
the atomic number of the constituent clement). 
This expression should be regarded as accurate 
to within a factor of a few while Eqs. (1) and 
especially T^^fj represent typical upper limits for 
the critical temperature. 

In Fig. 1 we have plotted the critical conden- 
sation temperatures for H and Fe, as a function 
of B, toghether with the coolest, thermally emit- 
ting NSs for which an estimate of the magnetic 
field (as computed from the spin-down formula) 
is available (see Tab. 1). We have also drawn 
an horizontal line at the surface temperature of 
RX J1856. 5-3754, for which B is unknown. We 
used as a color temperature the value derived 
from a blackbody fit in the X-rays, T^h- The sur- 
face temperature T,t„./ follows by applying the 
gravitational red-shift factor. 

As we can see from Fig. 1, most sources 
have Tsurf well in excess of the H critical tem- 
perature. Therefore, if surface layers are H- 
dominated, the presence of a gaseous atmo- 
sphere is ^inescapable. The only exception is 
RX J1308. 6-1-2127, for which an exceptionally 
strong field has been reported by [15]. How- 
ever, this value is probably preliminary and a re- 
analysis of the timing properties of this source is 
presently under way (Haberl, as reported at the 
34*'' COSPAR meeting). Similarly, H condensa- 
tion in RX J1856. 5-3754 requires the star to be a 
magnetar {B> lO" G). 

On the other hand, if NSs have not accreted 
much gas, we might detect thermal emission di- 
rectly from the iron surface layers. If this is 
the case, for B> 10^^ G the outermost layers of 
RX J1856. 5-3754 might be in the form of hot con- 
densed matter and the usual radiative transfer 
computations do not apply.^ The question of the 
nature of the emitted spectrum then arises. 

2.2. The Emitted Spectrum 

We consider a NS with metallic Fe surface lay- 
ers and repeat a computation, first carried out by 
[22] and based on the evaluation of the surface re- 
flectivity, for the parameter range relevant to cold 
isolated NSs. At each magnetic co-latitude 6 we 

^The case of RX ,10720.4-3125 is less certain, because the 
source falls only marginally within the region where Fe 
condensation is possible. 



Table 1 

Isolated Neutron Stars Parameters 



Source 


Tbb (eV) 


S(10i2 G) 


Refs." 


RX J1856.5-3754 


61.1 ±0.3 




[4,5] 


RX J0720.4-3125 


86.0 ±0.6 


2i.3i°:l 


[2,3] 


RX J1308.6-I-2127 


90.6 ±1.6 




[15] 


Vela 


128.4 ± 7 


3.3 


[16,17] 


Geminga 




1.5 


[18,19] 


PSR 0656-M4 


69.0 ±2.5 


4.7 


[16,20] 


PSR 1055-52 




1.1 


[16,21] 



flrst compute the total reflectivity pa,{i,/3,9) of 
the surface element for incident unpolarized radi- 
ation. Here i and P are the co-latitude and the az- 
imuthal angle of the incident wave vector, relative 
to the surface normal. The absorption coeflicient 
is then a^^ = 1 — p^, and Kirchhoff's law is used 
to derive the emissivity = ai^B^^{T), where T 
is the local temperature and B,^ is the blackbody 
function. The monochromatic flux dF^^ emitted 
by the surface element is computed by averaging 
over all incident directions, and further integra- 
tion over the star surface gives the total flux 

Fu, = / dF^= B^{T)sm6d6x 
Jo Jo 

/ / a^{i,P,e)smididl3. (2) 

Jo J-n/2 

In performing the last step we assumed a dipolar 
field, B ^ Sp[(4- /) cos2 6' + /]i/V2, where Bp is 
the polar field strength and / ~ 1.2 accounts for 
general-relativistic corrections in a Schwarzschild 
space-time. Wc computed two sets of models ei- 
ther by taking a constant surface temperature 
or allowing for a meridional temperature depen- 
dence as given by [23]. We also explored the pa- 
rameter space by varing the electron plasma fre- 
quency, uip, around = u}p{p = ps) , to account 
for possible deviations of the surface density from 
its zero-pressure value ps (see e.g. [13]). 

For B> 10^^ G, features related to the strong 
absorption around the electron cyclotron fre- 
quency fall well outside the X-ray range accessible 
to the Chandra LETGS and XMM-Newton EPIC- 
PN, and arc of no immediate interest. However, 
even above this field, we do not expect the spec- 
trum emitted by the solid to be, in general, a 
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Figure 1. Tcrit for H and Fc; as a function of B. 
Condensation is possible in the shaded region for 
Fe and in the cross-hatched region for H. Filled 
circles mark the position of the sources listed in 
Tab. 1; the horizontal line corresponds to the sur- 
face temperature of RX J1856. 5-3754. 



blackbody. For 100 eV, the mean energy of 
the refracted waves is comparable to (or lower 
than) the plasma frequency of the solid. There- 
fore, we encounter modes with very large re- 
fractivity (like whistlers in the terrestrial atmo- 
sphere), modes which are frozen in the medium, 
or modes highly damped which just enter the NS 
crust and disappear within a small penetration 
depth. All these effects may modify the spectrum 
through the production of features and edges, 
which in turn may or may not appear in the X- 
rays depending on the model parameters. 

In particular we found that, for i?p w 3 — 5 x 
10^^ G, the spectrum in the soft X-rays, (0.1- 
2 keV) band, shows negligible departures (within 
^ 10%) from the best-fitting blackbody and ex- 
hibits no features whatsoever. We point out that 
larger deviations are at the lowest energies, and 
that the agreement with the planckian is better 
than 2% if we restrict to the interval 0.14-2 keV. 
Spectra are shown in Fig. 2 for Bp = 3 x 10^^ G 
and Teff = 75 eV, which corresponds to 60 eV 
once the gravitational red-shift is applied. The 
constant temperature models have virtually no 
hardening, while we get T(,(,/Te// = 7 ~ 1.14 ac- 
counting for a meridional temperature variation. 
An important point is that, despite being close 
to planckian in shape, the spectrum is substan- 
tially depressed with respect to the blackbody at 



Figure 2. The spectrum at the star surface for 
Bp = 3 X 10^3 G, Teff = 75 eV and two val- 
ues of ojp-. Upfi (upper solid curve) and 2.5a;p,o 
(lower solid curve). Left and right panels refer 
to the uniform and meridional temperature dis- 
tributions, respectively. The dashed line is the 
blackbody at Teff and the dash-dotted line the 
best- fitting blackbody in the 0.1-2 keV range. 



Teff, owing to the reduced emissivity of the sur- 
face. The total power radiated by the star in the 
0.1-2 keV band is ~ 30-50% with respect to that 
of a perfect blackbody emitter and decreases with 
increasingly (Fig. 3). 

3. The Case of RX J1856.5-3754 

3.1. The X-ray Emission 

It has been recently proposed by a num- 
ber of authors ([10,7,11]) that the multiwave- 
length spectral energy distribution (SED) of 
RX J1856. 5-3754 may be explained in terms of 
a two-temperature surface distribution in a cool- 
ing NS. Although these models appear promis- 
ing, they are all based on the presumption that 
the emitted spectrum is planckian in shape. A 
number of mechanisms (magnetic smearing, ro- 
tation..) have been suggested in order to sup- 
press the spectral features predicted by atmo- 
spheric models, but no conclusive evidence has 
been provided yet that a nearly blackbody, fea- 
tureless spectrum can be emitted by an extended 
atmosphere covering the stellar crust. [11] discuss 
in detail the role of rotation, showing that phase- 
dependent Doppler shifts in a rapidly rotating 
neutron star (P w 1 ms) wash out all features, 
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Figure 3. Ratio of the emitted to the blackbody 
power in the 0.1-2 keV band for different values 
of the plasma frequency. Circles refer to Bp = 3 x 
1013 G and diamonds to Bp = 5 x lO^^ G. Filled 
and open symbols are for the uniform/meridional 
variation temperature distributions respectively. 



leaving a nearly planckian spectrum. Although 
such a short period can not be excluded on the ba- 
sis of present data, the detected periods of other 
thermally emitting INSs are in the range « 0.1- 
10 s, about two orders of magnitude larger. Also, 
in the picture by [11] the genuine surface tem- 
perature should correspond to the cooler black- 
body at T ~ 15 eV. Conventional cooling curves 
then imply a star's age of w 10® yr, which is 
hardly compatible with a ms period. Further- 
more, the energetics of the bow-shock nebula im- 

1 /2 

plies P = 4.6 G) ' ms. A ms spin period 

therefore necessary demands for a very low field 
star. Such a low field seems hard to reconcile with 
the limit on the age derived again from the bow 
shock energetics, (B/IO^^ G)(r/10'5 yr) ~ 3 - 4. 

Without a measured period and period deriva- 
tive, the magnetic field of RX J1856. 5-3754 is 
still a mystery. We found that, for the surface 
layers of RX J1856.5-3754 to be in the form of 
condensed Fe, its field should be in excess of 
10^3 G, and more probably at least 3-5 xlO^^ 
G (Fig. 1). Although rather high, this field 
strength is well below the magnetar range and 
is noticeably shared by another ROSAT isolated 
NS: RX J0720.4-3125 [3]. When combined with 
the limit derived from the bow shock energetics, 
(5/10^2 G)(r/106 yr) ~ 3 - 4, such field de- 
mands for a middle-aged NS (r » 10^ yr). The 



spectrum from a bare NS with a dipolar field 

in this range is shown in Fig. 2. In the 0.15- 
2 keV band the spectrum is featureless; devia- 
tions from the best-fitting blackbody are ^ 2%, 
well within the limits of calibration uncertainties 
of the Chandra LETGS (~ 10%, as reported by 
[11]). Our surface emissivity is angle-dependent, 
and we checked that no pulsations are expected 
to within 3%, the latest published value ([5]). 

For the radiation radius to be representative of 
the true star radius the emission properties of the 
surface need to be accounted for and, at least for 
thermal emission, we can use the spectral harden- 
ing 7 to quantify spectral deviations from a pure 
blackbody. Denoting the fraction of the stellar 
surface responsible for the observed emission and 
the ratio of emitted to blackbody power by /a 
and Je, the radiation radius is 



Rc 



:4.25/;-'/Vb'^V^^ioo 



Tbb 
60 eV 



km, (3) 



where dioo = d/(100pc). The "angular 
size" of RX J1856.5-3754 reported by [5] is 
i?oo/(rf/100pc) = 4.12 ± 0.68 km. If we correct 
this value by assuming d = 130 pc, emission from 
the entire surface {/a = 1) and uniform temper- 
ature (7 = 1), we get 7.7 ± 1.3 km< R^o^ 9.9 ± 
1.6 km in correspondence to 1 < LVp/iVpfi < 2.8 
(see Fig. 3). For the same range of plasma fre- 
quencies, models with meridional temperature 
variation give larger radii, between 9.4 ±1.5 km 
and 12 ± 2 km. Therefore, even this simple ap- 
proach can provide i?oc ^ 10-12 km, compatible 
with (soft) EOSs for M - 1.4Mq ([24]). 

3.2. The UV-Optical Excess 

The UV-optical flux from RX J1856.5-3754 is 

enhanced with respect to the extrapolation of the 
X-ray best-fitting blackbody. All proposed inter- 
pretations for the full SED (from optical to X- 
ray) require two components: either a cold sur- 
face plus an hotter cap or a bare quark star plus 
a thin crust with flux redistribution. Similarly, 
also our scenario requires the presence of a sec- 
ond component. A natural possibility is that the 
optical flux emitted from the NS surface is repro- 
cessed in a thin, ionized gas layer on the top of the 
Fe solid. Hydrogen can be deposited on the sur- 
face as the result of very slow accretion (w 10^ g 
in « 10'"' yr). The free- free absorption depth in 
the optical (~ 1-10 eV) is « 10*^ times larger 
than in the X-rays (^ 100 eV), thus in a range 
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of average densities (around « 10~^gcm~^) the 
gasoues layer is optically thin in the X-rays and 
thick to optical photons. The scale-height derived 
from hydrostatic equilibrium is « 0.1 — 1 cm, 
comparable with the scale-height at which elec- 
tron conduction from the crust keeps the gas al- 
most isothermal at the star surface temperature 
(« 10^ K). Preliminary calculation of the multi- 
wavelength SED are shown in Fig. 4. Parameters 
are the same as for the model with 2.5a;p,o shown 
in Fig. 3 and the density at the base of the layer 
is 2 X lO^'^gcm^'^. As expected the low-energy 
SED follows a Rayleigh-Jeans distribution at the 
gas (and star) temperature. With increasing en- 
ergy the optical depth of the layer drops below 
unity and therefore X-rays emitted by the solid 
surface emerge unchanged from the layer. How- 
ever, since the X-ray spectrum emitted by the 
star surface is depressed with respect to the black- 
body at the star temperature, the optical emis- 
sion appears enhanced. In the model presented 
here the optical- UV excess over the best-fitting X- 
ray blackbody is a factor ~ 3-4, somewhat lower 
than the value ~ 6 reported by [7] . Larger ratios 
might be obtained with a slightly larger gas tem- 
perature and a more thorough investigation is in 
progress. 

4. Conclusions 

We have shown that the observed X-ray-to- 
optical spectra of RX J1856. 5-3754 can be plau- 
sibly explained in terms of the emission from the 
surface of a "naked" neutron star. The absence 
of an atmosphere is due to a phase transition 
to a solid condensate, which can occur on cool 
(< 10^ K) neutron stars endowed with strong 
magnetic fields (J5> 10"'^^ G) and with metal- 
dominated outer layers. We have shown that, 
when B sa 10^^ G, the X-ray spectrum is fea- 
tureless and planckian in shape. The observed 
UV-optical enhanced emission can be explained 
by the presence of a gasoues, thin H shell, where 
the optical flux is reprocessed. 

While we caution that current limitations in 
our understanding of metallic condensates in 
strong magnetic fields renders our estimates of 
the surface reflectivity of bare NSs somewhat 
uncertain, we have shown that our model for 
RX J1856. 5-3754 predicts a value for the appar- 
ent radius of i?oo~ 12 km. For a canonical NS of 
mass 1.4 Mq, such a radius favors soft EOSs. 
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Figure 4. The multiwavclcngth SED of a bare 
neutron star covered by a thin H layer (full line); 
the dotted line is the spectrum emitted by the 
solid. Details as in Fig. 2. 
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